Overexpression of the UDP-glycosyltransferase UGT34A23 confers resistance to the diamide insecticide chlorantraniliprole in the tomato leafminer, Tuta absoluta.

The tomato leafminer, Tuta absoluta, is an invasive crop pest that has evolved resistance to many of the insecticides used for its control. To facilitate the investigation of the underpinning mechanisms of resistance in this species we generated a contiguous genome assembly using long-read sequencing data. We leveraged this genomic resource to investigate the genetic basis of resistance to the diamide insecticide chlorantraniliprole in Spanish strains of T. absoluta that exhibit high levels of resistance to this insecticide. Transcriptomic analyses revealed that, in these strains, resistance is not associated with previously reported target-site mutations in the diamide target-site, the ryanodine receptor, but rather is associated with the marked overexpression (20- to >100-fold) of a gene encoding a UDP-glycosyltransferase (UGT). Functional expression of this UGT, UGT34A23, via ectopic expression in Drosophila melanogaster demonstrated that it confers strong and significant resistance in vivo. The genomic resources generated in this study provide a powerful resource for further research on T. absoluta. Our findings on the mechanisms underpinning resistance to chlorantraniliprole will inform the development of sustainable management strategies for this important pest.

Differential regulation of the Tor gene homolog drives the red/green pigmentation phenotype in the aphid Myzuspersicae.

In some aphid species, intraspecific variation in body colour is caused by differential carotenoid content: whilst green aphids contain only yellow carotenoids (ß-, γ-, and ß,γ-carotenes), red aphids additionally possess red carotenoids (torulene and 3,4-didehydrolycopene). Unusually, within animals who typically obtain carotenoids from their diet, ancestral horizontal gene transfer of carotenoid biosynthetic genes from fungi (followed by gene duplication), have imbued aphids with the intrinsic gene repertoire necessary to biosynthesise carotenoids. In the pea aphid, Acyrthosiphon pisum a lycopene (phytoene) desaturase gene (Tor) underpins the red/green phenotype, with this locus present in heterozygous form in red individuals but absent in green aphids, resulting in them being unable to convert lycopene into the red compounds 3,4-didehydrolycopene and torulene. The green peach aphid, Myzus persicae, separated from the pea aphid for ≈45MY also exists as distinct colour variable morphs, with both red and green individuals present. Here, we examined genomic data for both red and green morphs of M. persicae and identified an enlarged (compared to A. pisum) repertoire of 16 carotenoid biosynthetic genes (11 carotenoid desaturases and five carotenoid cyclase/synthase genes). From these, we identify the homolog of A. pisum Tor (here called carotene desaturase 2 or CDE-2) and show through 3D modelling that this homolog can accommodate the torulene precursor lycopene and, through RNA knockdown feeding experiments, demonstrate that disabling CDE-2 expression in red M. persicae clones results in green-coloured offspring. Unlike in A. pisum, we show that functional CDE-2 is present in the genomes of both red and green aphids. However, expression differences between the two colour morphs (350-700 fold CDE-2 overexpression in red clones), potentially driven by variants identified in upstream putative regulatory elements, underpin this phenotype. Thus, whilst aphids have a common origin of their carotenoid biosynthetic pathway, two aphid species separated for over 40MY have evolved very different drivers of intraspecific colour variation.

The genome sequence of the Neotropical brown stink bug, Euschistus heros provides insights into population structure, demographic history and signatures of adaptation.

The Neotropical brown stink bug, Euschistus heros, is a major pest of soybean in South America. The importance of E. heros as a pest has grown significantly in recent times due to increases in its abundance and range, and the evolution of insecticide resistance. Recent work has begun to examine the genetic diversity, population structure, and genetic mechanisms of insecticide resistance in E. heros. However, to date, investigation of these topics has been hampered by a lack of genomic resources for this species. Here we address this need by assembling a high-quality draft genome for E. heros. We used a combination of short and long read sequencing to assemble an E. heros genome of 1.4 Gb comprising 906 contigs with a contig N50 of 3.5 MB. We leveraged this new genomic resource, in combination with genotyping by sequencing, to explore genetic diversity in populations of this species in Brazil and identify genetic loci in the genome which are under selection. Our genome-wide analyses, confirm that there are two populations of E. heros co-occurring in different geographical regions in Brazil, and that, in certain regions of the country these populations are hybridizing. We identify several regions of the genome as under selection, including markers associated with putative insecticide resistance genes. Taken together, the new genomic resources generated in this study will accelerate research into fundamental aspects of stinkbug biology and applied aspects relating to the sustainable control of a highly damaging crop pest.

Integrative omics approaches for biosynthetic pathway discovery in plants.

Covering: up to 2022With the emergence of large amounts of omics data, computational approaches for the identification of plant natural product biosynthetic pathways and their genetic regulation have become increasingly important. While genomes provide clues regarding functional associations between genes based on gene clustering, metabolome mining provides a foundational technology to chart natural product structural diversity in plants, and transcriptomics has been successfully used to identify new members of their biosynthetic pathways based on coexpression. Thus far, most approaches utilizing transcriptomics and metabolomics have been targeted towards specific pathways and use one type of omics data at a time. Recent technological advances now provide new opportunities for integration of multiple omics types and untargeted pathway discovery. Here, we review advances in plant biosynthetic pathway discovery using genomics, transcriptomics, and metabolomics, as well as recent efforts towards omics integration. We highlight how transcriptomics and metabolomics provide complementary information to link genes to metabolites, by associating temporal and spatial gene expression levels with metabolite abundance levels across samples, and by matching mass-spectral features to enzyme families. Furthermore, we suggest that elucidation of gene regulatory networks using time-series data may prove useful for efforts to unwire the complexities of biosynthetic pathway components based on regulatory interactions and events.

Genome assembly of Danaus chrysippus and comparison with the Monarch Danaus plexippus.

Milkweed butterflies in the genus Danaus are studied in a diverse range of research fields including the neurobiology of migration, biochemistry of plant detoxification, host-parasite interactions, evolution of sex chromosomes, and speciation. We have assembled a nearly chromosomal genome for Danaus chrysippus (known as the African Monarch, African Queen, and Plain Tiger) using long-read sequencing data. This species is of particular interest for the study of genome structural change and its consequences for evolution. Comparison with the genome of the North American Monarch Danaus plexippus reveals generally strong synteny but highlights 3 inversion differences. The 3 chromosomes involved were previously found to carry peaks of intraspecific differentiation in D. chrysippus in Africa, suggesting that these inversions may be polymorphic and associated with local adaptation. The D. chrysippus genome is over 40% larger than that of D. plexippus, and nearly all of the additional ∼100 Megabases of DNA comprises repeats. Future comparative genomic studies within this genus will shed light on the evolution of genome architecture.

Global patterns in genomic diversity underpinning the evolution of insecticide resistance in the aphid crop pest Myzus persicae.

The aphid Myzus persicae is a destructive agricultural pest that displays an exceptional ability to develop resistance to both natural and synthetic insecticides. To investigate the evolution of resistance in this species we generated a chromosome-scale genome assembly and living panel of >110 fully sequenced globally sampled clonal lines. Our analyses reveal a remarkable diversity of resistance mutations segregating in global populations of M. persicae. We show that the emergence and spread of these mechanisms is influenced by host-plant associations, uncovering the widespread co-option of a host-plant adaptation that also offers resistance against synthetic insecticides. We identify both the repeated evolution of independent resistance mutations at the same locus, and multiple instances of the evolution of novel resistance mechanisms against key insecticides. Our findings provide fundamental insights into the genomic responses of global insect populations to strong selective forces, and hold practical relevance for the control of pests and parasites.

Molecular innovations underlying resistance to nicotine and neonicotinoids in the aphid Myzus persicae.

The green peach aphid, Myzus persicae, is a globally distributed highly damaging crop pest. This species has demonstrated an exceptional ability to evolve resistance to both synthetic insecticides used for control, and natural insecticides produced by certain plants as a chemical defense against insect attack. Here we review work characterizing the evolution of resistance in M. persicae to the natural insecticide nicotine and the structurally related class of synthetic neonicotinoid insecticides. We outline how research on this topic has provided insights into long-standing questions of both evolutionary and applied importance. These include questions pertaining to the origins of novel traits, the number and nature of mutational events or 'adaptive steps' underlying the evolution of new phenotypes, and whether host plant adaptations can be co-opted to confer resistance to synthetic insecticides. Finally, research on the molecular mechanisms underlying insecticide resistance in M. persicae has generated several outstanding questions on the genetic architecture of resistance to both natural and synthetic xenobiotics, and we conclude by identifying key knowledge gaps for future research. © 2021 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

The genetic architecture of a host shift: An adaptive walk protected an aphid and its endosymbiont from plant chemical defenses.

Host shifts can lead to ecological speciation and the emergence of new pests and pathogens. However, the mutational events that facilitate the exploitation of novel hosts are poorly understood. Here, we characterize an adaptive walk underpinning the host shift of the aphid Myzus persicae to tobacco, including evolution of mechanisms that overcame tobacco chemical defenses. A series of mutational events added as many as 1.5 million nucleotides to the genome of the tobacco-adapted subspecies, M. p. nicotianae, and yielded profound increases in expression of an enzyme that efficiently detoxifies nicotine, both in aphid gut tissue and in the bacteriocytes housing the obligate aphid symbiont Buchnera aphidicola. This dual evolutionary solution overcame the challenge of preserving fitness of a mutualistic symbiosis during adaptation to a toxic novel host. Our results reveal the intricate processes by which genetic novelty can arise and drive the evolution of key innovations required for ecological adaptation.

De Novo Genome Assembly of the Meadow Brown Butterfly, Maniola jurtina.

Meadow brown butterflies (Maniola jurtina) on the Isles of Scilly represent an ideal model in which to dissect the links between genotype, phenotype and long-term patterns of selection in the wild - a largely unfulfilled but fundamental aim of modern biology. To meet this aim, a clear description of genotype is required. Here we present the draft genome sequence of M. jurtina to serve as a founding genetic resource for this species. Seven libraries were constructed using pooled DNA from five wild caught spotted females and sequenced using Illumina, PacBio RSII and MinION technology. A novel hybrid assembly approach was employed to generate a final assembly with an N50 of 214 kb (longest scaffold 2.9 Mb). The sequence assembly described here predicts a gene count of 36,294 and includes variants and gene duplicates from five genotypes. Core BUSCO (Benchmarking Universal Single-Copy Orthologs) gene sets of Arthropoda and Insecta recovered 90.5% and 88.7% complete and single-copy genes respectively. Comparisons with 17 other Lepidopteran species placed 86.5% of the assembled genes in orthogroups. Our results provide the first high-quality draft genome and annotation of the butterfly M. jurtina.

Whole-chromosome hitchhiking driven by a male-killing endosymbiont.

Neo-sex chromosomes are found in many taxa, but the forces driving their emergence and spread are poorly understood. The female-specific neo-W chromosome of the African monarch (or queen) butterfly Danaus chrysippus presents an intriguing case study because it is restricted to a single 'contact zone' population, involves a putative colour patterning supergene, and co-occurs with infection by the male-killing endosymbiont Spiroplasma. We investigated the origin and evolution of this system using whole genome sequencing. We first identify the 'BC supergene', a broad region of suppressed recombination across nearly half a chromosome, which links two colour patterning loci. Association analysis suggests that the genes yellow and arrow in this region control the forewing colour pattern differences between D. chrysippus subspecies. We then show that the same chromosome has recently formed a neo-W that has spread through the contact zone within approximately 2,200 years. We also assembled the genome of the male-killing Spiroplasma, and find that it shows perfect genealogical congruence with the neo-W, suggesting that the neo-W has hitchhiked to high frequency as the male-killer has spread through the population. The complete absence of female crossing-over in the Lepidoptera causes whole-chromosome hitchhiking of a single neo-W haplotype, carrying a single allele of the BC supergene and dragging multiple non-synonymous mutations to high frequency. This has created a population of infected females that all carry the same recessive colour patterning allele, making the phenotypes of each successive generation highly dependent on uninfected male immigrants. Our findings show how hitchhiking can occur between the physically unlinked genomes of host and endosymbiont, with dramatic consequences.

Host plant adaptation in the polyphagous whitefly, Trialeurodes vaporariorum, is associated with transcriptional plasticity and altered sensitivity to insecticides.

BACKGROUND: The glasshouse whitefly, Trialeurodes vaporariorum, is a damaging crop pest and an invasive generalist capable of feeding on a broad range of host plants. As such this species has evolved mechanisms to circumvent the wide spectrum of anti-herbivore allelochemicals produced by its host range. T. vaporariorum has also demonstrated a remarkable ability to evolve resistance to many of the synthetic insecticides used for control. RESULTS: To gain insight into the molecular mechanisms that underpin the polyphagy of T. vaporariorum and its resistance to natural and synthetic xenobiotics, we sequenced and assembled a reference genome for this species. Curation of genes putatively involved in the detoxification of natural and synthetic xenobiotics revealed a marked reduction in specific gene families between this species and another generalist whitefly, Bemisia tabaci. Transcriptome profiling of T. vaporariorum upon transfer to a range of different host plants revealed profound differences in the transcriptional response to more or less challenging hosts. Large scale changes in gene expression (> 20% of genes) were observed during adaptation to challenging hosts with a range of genes involved in gene regulation, signalling, and detoxification differentially expressed. Remarkably, these changes in gene expression were associated with significant shifts in the tolerance of host-adapted T. vaporariorum lines to natural and synthetic insecticides. CONCLUSIONS: Our findings provide further insights into the ability of polyphagous insects to extensively reprogram gene expression during host adaptation and illustrate the potential implications of this on their sensitivity to synthetic insecticides.

The leafcutter bee, Megachile rotundata, is more sensitive to N-cyanoamidine neonicotinoid and butenolide insecticides than other managed bees.

Recent research has shown that several managed bee species have specific P450 enzymes that are preadapted to confer intrinsic tolerance to some insecticides including certain neonicotinoids. However, the universality of this finding across managed bee pollinators is unclear. Here we show that the alfalfa leafcutter bee, Megachile rotundata, lacks such P450 enzymes and is >2,500-fold more sensitive to the neonicotinoid thiacloprid and 170-fold more sensitive to the butenolide insecticide flupyradifurone than other managed bee pollinators. These findings have important implications for the safe use of insecticides in crops where M. rotundata is used for pollination, and ensuring that regulatory pesticide risk assessment frameworks are protective of this species.

Phylogenomic analyses and distribution of terpene synthases among Streptomyces.

Terpene synthases are widely distributed among microorganisms and have been mainly studied in members of the genus Streptomyces. However, little is known about the distribution and evolution of the genes for terpene synthases. Here, we performed whole-genome based phylogenetic analysis of Streptomyces species, and compared the distribution of terpene synthase genes among them. Overall, our study revealed that ten major types of terpene synthases are present within the genus Streptomyces, namely those for geosmin, 2-methylisoborneol, epi-isozizaene, 7-epi-α-eudesmol, epi-cubenol, caryolan-1-ol, cyclooctat-9-en-7-ol, isoafricanol, pentalenene and α-amorphene. The Streptomyces species divide in three phylogenetic groups based on their whole genomes for which the distribution of the ten terpene synthases was analysed. Geosmin synthases were the most widely distributed and were found to be evolutionary positively selected. Other terpene synthases were found to be specific for one of the three clades or a subclade within the genus Streptomyces. A phylogenetic analysis of the most widely distributed classes of Streptomyces terpene synthases in comparison to the phylogenomic analysis of this genus is discussed.

Genomic insights into neonicotinoid sensitivity in the solitary bee Osmia bicornis.

The impact of pesticides on the health of bee pollinators is determined in part by the capacity of bee detoxification systems to convert these compounds to less toxic forms. For example, recent work has shown that cytochrome P450s of the CYP9Q subfamily are critically important in defining the sensitivity of honey bees and bumblebees to pesticides, including neonicotinoid insecticides. However, it is currently unclear if solitary bees have functional equivalents of these enzymes with potentially serious implications in relation to their capacity to metabolise certain insecticides. To address this question, we sequenced the genome of the red mason bee, Osmia bicornis, the most abundant and economically important solitary bee species in Central Europe. We show that O. bicornis lacks the CYP9Q subfamily of P450s but, despite this, exhibits low acute toxicity to the N-cyanoamidine neonicotinoid thiacloprid. Functional studies revealed that variation in the sensitivity of O. bicornis to N-cyanoamidine and N-nitroguanidine neonicotinoids does not reside in differences in their affinity for the nicotinic acetylcholine receptor or speed of cuticular penetration. Rather, a P450 within the CYP9BU subfamily, with recent shared ancestry to the Apidae CYP9Q subfamily, metabolises thiacloprid in vitro and confers tolerance in vivo. Our data reveal conserved detoxification pathways in model solitary and eusocial bees despite key differences in the evolution of specific pesticide-metabolising enzymes in the two species groups. The discovery that P450 enzymes of solitary bees can act as metabolic defence systems against certain pesticides can be leveraged to avoid negative pesticide impacts on these important pollinators.

Extension of Partial Gene Transcripts by Iterative Mapping of RNA-Seq Raw Reads.

Many non-model organisms lack reference genomes and the sequencing and de novo assembly of an organism's transcriptome is an affordable means by which to characterize the coding component of its genome. Despite the advances that have made this possible, assembling a transcriptome without a known reference usually results in a collection of full-length and partial gene transcripts. The downstream analysis of genes represented as partial transcripts then often requires further experimental work in the laboratory in order to obtain full- length sequences. We have explored whether partial transcripts, encoding genes of interest present in de novo assembled transcriptomes of a model and non-model insect species, could be further extended by iterative mapping against the raw transcriptome sequencing reads. Partial sequences encoding cytochrome P450s and carboxyl/cholinesterase were used in this analysis because they are large multigene families and exhibit significant variation in expression. We present an effective method to improve the continuity of partial transcripts in silico that, in the absence of a reference genome, maybe a quick and cost-effective alternative to their extension by laboratory experimentation. Our approach resulted in the successful extension of incompletely assembled transcripts, often to full length. We experimentally validated these results \textit{in silico} and using real-time PCR and sequencing.

Unravelling the Molecular Determinants of Bee Sensitivity to Neonicotinoid Insecticides.

The impact of neonicotinoid insecticides on the health of bee pollinators is a topic of intensive research and considerable current debate [1]. As insecticides, certain neonicotinoids, i.e., N-nitroguanidine compounds such as imidacloprid and thiamethoxam, are as intrinsically toxic to bees as to the insect pests they target. However, this is not the case for all neonicotinoids, with honeybees orders of magnitude less sensitive to N-cyanoamidine compounds such as thiacloprid [2]. Although previous work has suggested that this is due to rapid metabolism of these compounds [2-5], the specific gene(s) or enzyme(s) involved remain unknown. Here, we show that the sensitivity of the two most economically important bee species to neonicotinoids is determined by cytochrome P450s of the CYP9Q subfamily. Radioligand binding and inhibitor assays showed that variation in honeybee sensitivity to N-nitroguanidine and N-cyanoamidine neonicotinoids does not reside in differences in their affinity for the receptor but rather in divergent metabolism by P450s. Functional expression of the entire CYP3 clade of P450s from honeybees identified a single P450, CYP9Q3, that metabolizes thiacloprid with high efficiency but has little activity against imidacloprid. We demonstrate that bumble bees also exhibit profound differences in their sensitivity to different neonicotinoids, and we identify CYP9Q4 as a functional ortholog of honeybee CYP9Q3 and a key metabolic determinant of neonicotinoid sensitivity in this species. Our results demonstrate that bee pollinators are equipped with biochemical defense systems that define their sensitivity to insecticides and this knowledge can be leveraged to safeguard bee health.

Neofunctionalization of Duplicated P450 Genes Drives the Evolution of Insecticide Resistance in the Brown Planthopper.

Gene duplication is a major source of genetic variation that has been shown to underpin the evolution of a wide range of adaptive traits [1, 2]. For example, duplication or amplification of genes encoding detoxification enzymes has been shown to play an important role in the evolution of insecticide resistance [3-5]. In this context, gene duplication performs an adaptive function as a result of its effects on gene dosage and not as a source of functional novelty [3, 6-8]. Here, we show that duplication and neofunctionalization of a cytochrome P450, CYP6ER1, led to the evolution of insecticide resistance in the brown planthopper. Considerable genetic variation was observed in the coding sequence of CYP6ER1 in populations of brown planthopper collected from across Asia, but just two sequence variants are highly overexpressed in resistant strains and metabolize imidacloprid. Both variants are characterized by profound amino-acid alterations in substrate recognition sites, and the introduction of these mutations into a susceptible P450 sequence is sufficient to confer resistance. CYP6ER1 is duplicated in resistant strains with individuals carrying paralogs with and without the gain-of-function mutations. Despite numerical parity in the genome, the susceptible and mutant copies exhibit marked asymmetry in their expression with the resistant paralogs overexpressed. In the primary resistance-conferring CYP6ER1 variant, this results from an extended region of novel sequence upstream of the gene that provides enhanced expression. Our findings illustrate the versatility of gene duplication in providing opportunities for functional and regulatory innovation during the evolution of an adaptive trait.

SaDA: From Sampling to Data Analysis-An Extensible Open Source Infrastructure for Rapid, Robust and Automated Management and Analysis of Modern Ecological High-Throughput Microarray Data.

One of the most crucial characteristics of day-to-day laboratory information management is the collection, storage and retrieval of information about research subjects and environmental or biomedical samples. An efficient link between sample data and experimental results is absolutely important for the successful outcome of a collaborative project. Currently available software solutions are largely limited to large scale, expensive commercial Laboratory Information Management Systems (LIMS). Acquiring such LIMS indeed can bring laboratory information management to a higher level, but most of the times this requires a sufficient investment of money, time and technical efforts. There is a clear need for a light weighted open source system which can easily be managed on local servers and handled by individual researchers. Here we present a software named SaDA for storing, retrieving and analyzing data originated from microorganism monitoring experiments. SaDA is fully integrated in the management of environmental samples, oligonucleotide sequences, microarray data and the subsequent downstream analysis procedures. It is simple and generic software, and can be extended and customized for various environmental and biomedical studies.

Molecular tools for the selective detection of nine diatom species biomarkers of various water quality levels.

Our understanding of the composition of diatom communities and their response to environmental changes is currently limited by laborious taxonomic identification procedures. Advances in molecular technologies are expected to contribute more efficient, robust and sensitive tools for the detection of these ecologically relevant microorganisms. There is a need to explore and test phylogenetic markers as an alternative to the use of rRNA genes, whose limited sequence divergence does not allow the accurate discrimination of diatoms at the species level. In this work, nine diatom species belonging to eight genera, isolated from epylithic environmental samples collected in central Italy, were chosen to implement a panel of diatoms covering the full range of ecological status of freshwaters. The procedure described in this work relies on the PCR amplification of specific regions in two conserved diatom genes, elongation factor 1-a (eEF1-a) and silicic acid transporter (SIT), as a first step to narrow down the complexity of the targets, followed by microarray hybridization experiments. Oligonucleotide probes with the potential to discriminate closely related species were designed taking into account the genetic polymorphisms found in target genes. These probes were tested, refined and validated on a small-scale prototype DNA chip. Overall, we obtained 17 highly specific probes targeting eEF1-a and SIT, along with 19 probes having lower discriminatory power recognizing at the same time two or three species. This basic array was validated in a laboratory setting and is ready for tests with crude environmental samples eventually to be scaled-up to include a larger panel of diatoms. Its possible use for the simultaneous detection of diatoms selected from the classes of water quality identified by the European Water Framework Directive is discussed.

Molecular detection of a potentially toxic diatom species.

A few diatom species produce toxins that affect human and animal health. Among these, members of the Pseudo-nitzschia genus were the first diatoms unambiguously identified as producer of domoic acid, a neurotoxin affecting molluscan shell-fish, birds, marine mammals, and humans. Evidence exists indicating the involvement of another diatom genus, Amphora, as a potential producer of domoic acid. We present a strategy for the detection of the diatom species Amphora coffeaeformis based on the development of species-specific oligonucleotide probes and their application in microarray hybridization experiments. This approach is based on the use of two marker genes highly conserved in all diatoms, but endowed with sufficient genetic divergence to discriminate diatoms at the species level. A region of approximately 450 bp of these previously unexplored marker genes, coding for elongation factor 1-a (eEF1-a) and silicic acid transporter (SIT), was used to design oligonucleotide probes that were tested for specificity in combination with the corresponding fluorescently labeled DNA targets. The results presented in this work suggest a possible use of this DNA chip technology for the selective detection of A. coffeaeformis in environmental settings where the presence of this potential toxin producer may represent a threat to human and animal health. In addition, the same basic approach can be adapted to a wider range of diatoms for the simultaneous detection of microorganisms used as biomarkers of different water quality levels.